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Abstract 

We solve the operator ordering problem for the quantum continuous integrable su(l,l) 
Landau-Lifshitz model, and give a prescription to obtain the quantum trace identities, and 
the spectrum for the higher-order local charges. We also show that this method, based on op- 
erator regularization and renormalization, which guarantees quantum integrability, as well as 
the construction of self-adjoint extensions, can be used as an alternative to the discretization 
procedure, and unlike the latter, is based only on integrable representations. 
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1 Introduction 

It is well-known that the continuous integrable models are much harder to quantize than the 
discrete ones. On the other hand, the discretization of a given continuous model may not be 
known, and the standard techniques to obtain one lead to quite complicated results [THS]. One 
immediately faces several problems working directly with a continuous integrable model. The main 
difficulty is the problem of singularities in the operator product in the quantum theory, and, as a 
result, the operator ordering problem. There are several known continuous models that work out 
without any serious problems, e.g., the non-linear Schrodinger and Thirring models [4]. However, 
as it was emphasized in the previous work [7|, the derivation of the integrability in those models 
is not mathematically rigorous, although the formal manipulations with singular objects lead to 
the correct results. The mathematically strict method to deal with such models was first obtained 
in [3, where on the example of the highly singular sn(l, 1) Landau-Lifshitz (LL) model it was 
shown that not only one cannot do formal manipulations with singular potentials of the 6"{x) type, 
but it is necessary to construct self-adjoint extensions and work with regularized operator product. 
The physical consequences of such operations are significant. For example, the thermodynamical 
properties of the model considered in [8] turn out to be quite non-trivial, and allowed us to show 
the stability of the system. Another more important consequence concerns the widely accepted 
representation of the n-particle wave-function in the "factorizable" form, e.g., for n = 2: 



which, as we showed, is not correct in general. It is certainly not correct for the Landau-Lifshitz 
model, or the fermionic AAF model [^. It is also not correct, for the same reason, for the principal 
chiral model, quantization and integrability of which was considered in |10H12| . even though, in 
this case one does not encounter singular potentials such as 6"{x). Nevertheless, even in this usual 
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6{x) interaction it is necessary to construct self-adjoint extensions. Thus, the continuous integrable 
models should be treated with more mathematical care0 

In [7] we have worked out this method in details on the example of the quantum Hamiltonian for 
the su{l, 1) Landau-Lifshitz model. The natural question is how to extend this method to higher 
order charges. In this article we address this issue, and outline the procedure to obtain the quantum 
trace identities, which should be applicable for a wide range of continuous integrable models. Let 
us remind, that, for example, even for the much simpler NLS model, obtaining the quantum trace 
identities is not a simple task [TH5l [T3l^l5| . In particular, the widely accepted notion that the 
normal ordering with respect to the bosonic fields 'iIj{x) and ■ip~^{x) solves the operator ordering 
problem, is not quite correct. In fact, although there are no problems for the first three charges 
following from the quantum trace identities, the forth charge does not correspond to the usual 
normal ordering [ll[T6l[T7]. For the Landau-Lifshitz model, the situation is more complex, because 
now there are three fields S subject to be constrained on a hyperboloid. It is not clear a priori how 
to choose consistent operator ordering in this case. In the preliminary paper of Sklyanin |18| . the 
attempt was made to resolve the constraint, and then use the normal ordering for the remaining two 
operators. However, this procedure faces significant difficulties on the later stage. For example, in 
this formalism, it was not possible to obtain the quantum Hamiltonian in terms of the unconstrained 
fields, while in terms of the S fields the Hamiltonian, as we had shown in [7] has a quite simple 
form. 

This issue becomes even more important for other, more complicated continuous models, where 
resolving the constraints may not be technically possible. Thus, the first problem we address in 
this paper is how to work with operator product in this case. This is done with a special operator 
regularization and renormalization procedure. Namely, we show that it leads to the intertwining 
relation, which is fundamental for quantum integrability and, therefore, quantum trace identities. 
The regularization and renormalization of operators that we use essentially differ from the one 
used initially by Sklyanin. The latter was used to only satisfy the intertwining relation, thus, 
regularizing the singularities only of a special type. However, as we show below, it is not enough to 
correctly quantize the model, and, therefore, we introduce a more general operator regularization 
and renormalization, which takes care of all singularities, and which has a number of interesting 
properties and consequences. In particular, the regularized and renormalized operator algebra turns 
out to be highly non-local even for the isotropic LL-modelH and as a result, the subsequent analysis 
becomes much more involved. Here we also note, that the standard derivation of the quantum trace 
identities for continuous models does in fact use [4tll7|. as an intermediate step, the discretization of 
the model, which, as we emphasized above, is not always an obvious procedure. We discuss the steps 
to avoid this difficulty and at no point introduce any discretization. Moreover, we show that unlike 
the discretization procedure, which should utilize non-integrable representations, in our method we 
use only integrable representations, and verify this result for the next higher-order charge. 

The second step we consider is the necessary self-adjointness check of each operator, associ- 
ated with the conserved quantity, which requires the construction of self-adjoint extensions. This 

^ These different integrable models arise in different sub-sectors of the AdSs x string and correspondingly 
require different constructions of self-adjoint extensions. The interesting question of constructing a larger space to 
accommodate the full string sector will be investigated in a separate publication. 

^ Let us note, that such non-local algebras (the Sklyanin algebras) originally arose from the anisotropic LL 
model |18ti20| . Here we show that it is necessary to introduce an analogue of the Sklyanin algebra even for the 
isotropic case. 
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leads to mathematically correct description in the quantum mechanical formalism [7], and which 
clearly shows that the intuitively obvious expressions of the type (jl.ip are not always correct. We 
demonstrate this check for the Hamiltonian, and the next non-trivial cubic operator. 

The third step in our method utilizes the A^-operator formalism, which establishes a connection 
between our method and the quantum inverse scattering method, and which makes it possible 
to actually derive the spectrum for each operator without, as it was initially done by Sklyanin, 
guessing it and using its classical limit. We again demonstrate this point for the Hamiltonian and 
the next non-trivial charge. In the process, we also discuss a number of important related points. 
In particular, the constraint regularization and renormalization are carefully considered, and the 
quantum analogue of it is derived. We show that after all regularizations are removed, the resulting 
Casimir operator corresponds to integrable representations, consistent with self-adjointness, and 
boundary conditions. 

We also present the classical inverse scattering method analysis for the sn(l, 1) LL-model, on 
which we subsequently rely during the analysis of the quantum model. Finally, by comparing our 
results to the classical ones, and using the symmetry of the model, we argue that any higher- 
order quantum conserved charge can be written in terms of the three fundamental charges: the 
Casimir operator, the quantum Hamiltonian, and the next non-trivial conserved local charge, which 
essentially solves the problem of quantum trace identities for the su(l, 1) continuous LL-model. 

The paper is organized as follows. In Section [2] we give the main relevant expressions and 
definitions related to the quantization of the su{l, 1) LL-model. In particular, in Section [2. II we give 
all essential results of the classical inverse scattering method. In Sections 12.21 and 12.31 we overview 
the quantization procedure, and discuss the Sklyanin product, as well as the construction of self- 
adjoint extensions and the hierarchy of A^m-operators. The main results are presented in Section [3] 
where we introduce the operator regularization, discuss the resulting algebraic structure and show 
the quantum integrability. We also derive the quantum fundamental charges, their spectrum, and 
establish a connection to the quantum inverse scattering method. Finally, in Section |4] we check 
the self-adjointness of the higher-order fundamental charge. In conclusion. Section [5l we discuss 
the obtained results, and propose further directions to develop the ideas presented in this paper for 
other, more interesting continuous integrable models. 

2 Overview of the Landau-Lifshitz model 

In this section we make a short review, following |18 |[2Tl - [23] . of the Landau-Lifshitz model in 
the context of the inverse scattering method and discuss the arising difficulties upon the attempt 
to obtain algebraic expressions for the conserved charges in the quantum theory. 

The most general anisotropic Landau-Lifshitz model is described in the classical case by the 
following Hamiltonian (see Section [2. II for more details): 





where e = ±1 and the fields {i = 1,2,3) define a vector in three-dimensional space, S 
(5"^, (S*^, 5^) with the scalar product defined by: 



(2.2) 
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It is convenient to introduce = i: iS'^, in terms of which the Poisson structure has the form: 

{S\x),S^{y)} = ±iS^ix)5{x-y) 

(2.3) 

{S'ix),S^iy)} = 2ieS''ix)5ix-y) 

The sign of £ determines whether the model is defined on the su{2) sphere (e = —1) or on the 
1) hyperboloid (e = 1). As noted in |18| . it is only possible to construct physically meaningful 
states over the ferromagnetic vacuum for the hyperboloid case. The role of the parameter 7 is 
to introduce anisotropy in the model and, as explained in |18| . the isotropic case (7 = 0) should 
be treated separately from the more general anisotropic (7 7^ 0) one. Unlike the former case, 
where the standard inverse scattering method works out throughout, that is, the Yang-Baxter and 
bilinear equations are satisfied with the usual choice of the i2-matrix and yield the commuting family 
of operators, the latter requires a modification of the monodromy matrix and the spin operator 
algebra (giving rise to the quadratic Sklyanin algebra |19|I20|). so that the intertwining relations 
have a solution with the i2-matrix of the XyZ-model. In this article we focus only on the isotropic 
su{l, 1) case. 



2.1 Classical inverse scattering method 

We start our brief review of the classical inverse scattering method with fixing the notations. 
This will allow us to treat both su{2) and su{l, 1) models simultaneously (see below the crucial 
difference between these models in the quantum case). Introducing the metric 

rjij = diag{-£, -e, 1) (2.4) 

where again e = —1 corresponds to the su{2) case and e = 1 corresponds to the su{l, 1) one, we 
can write the equations of motion for the spin field, S^{x), in the usual form: 

d S'* 

— = e'^'^Sjd^Sk (2.5) 

where e^"^^ = 1, d := and we used the metric (12. 4p to lower indices, i.e.. Si = rjijS^ . It is clear 
that the dynamics is consistent with the constraint 

VijS'S^ = 1 (2.6) 

The equation (j2.5p plus the constraint (I2.6P define the classical dynamics of the su{2) and su(l, 1) 
models. 

It is easily verified that (|2.5p can be obtained as the canonical equations using the following 
Hamiltonian and Poisson structure (c.f., Eqs. (|2.m2.3p ): 

H = -^jdx riijdS\x)dS^ {x) (2.7) 
{S\x),S^{y)}=ee'^''Sk{y)6{x-y) (2.8) 
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Now we present a brief account on how to solve (|2.5p using the inverse scattering method. Here we 
mostly follow |21fl23) . generalizing it to incorporate su(l, 1) case. With the help of cr-matrices 

[a\a^] = -2ie€'^''ak (2.9) 
{a\a^} = 2ri'^ (2.10) 

one can introduce a matrix- valued field S: 

S := 5Vi (2.11) 
It can be easily checked that Eqs. (|2.5|) and (|2.6|) can be re-written as follows: 

^=4[S,#S1 (2.2) 

= 1 (2.13) 

To pass to the auxiliary spectral problem, we need to construct a Lax pair. By a direct calculation, 
it is verified that the following choice 

L = Sd and M = -e{2Sd'^ + dSd) (2.14) 



makes (|2.5|) and (|2.6|) equivalent to: 

i) the auxiliary spectral problem: 

ii) the time evolution of ip, given by: 



iLip = Xi/j (2.15) 



i^=Mip (2.16) 

As a result, the crucial property, that A is a constant of motion, can be proven. 

With the help of the constraint (j2.6p the spectral equation (|2.15p can be brought to a more 
convenient forn^fl 

idi; = XS^j (2.17) 
It is important to notice that this equation as well as the asymptotic for S: S(a;) — > do not 

\x\—>oo 

depend on e, i.e., both sets of invariants, for su{2) and su(l, 1), will be exactly the same if written 
in a covariant form, i.e., using the metric (12. 4p to contract indices (see below). 

Equations (j2.16p and (|2.17p are the starting point of the classical inverse scattering method. 
Here we only sketch the main steps leading to the construction of the infinite set of invariants. 

One starts with (|2.17p to construct two Jost solutions, F{x,X) and G{x,X), which satisfy the 
following asymptotical conditions: 

F{x,X) e"*^"'^ (2.18) 



x—^+oo 



G(x,A) e"*^"'^ (2.19) 

z— >--oo 



■^Thus, the classical L'^'^^^-operator, which enters the equation = L^'^'^^ip, has the form: 

S'-^{x) iS-{x) 



iS+{x) -S^{x) 
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These solutions depend on S^{x) through 8(2;) which plays the role of the potential in (j2.17p . 
Then the solution for the spin- field, 5*(x), is obtained by solving Gelfand-Marchenko equation, 
which though being integral equation, still is simpler than the original Landau-Lifshitz non-linear 
problem. 

The time dependence of the solutions is controlled by (j2.16p . It is important to notice that 
e in the definition of M-operator can be eliminated by the time redefinition: t — )• et. Then the 
evolution of both, su{2) and sn(l, 1), models is the same in terms of the new time. In particular, 
the invariants will be the same. Introducing the transition matrix, T{\), which connects two Jost 
solutions 

G{x,\) =F{x,X)T{X) 
and using (j2.16p one can easily find the time dependence of T(A): 

r(A, t) = e-2*^^''^'*r(A, 0)6^^^^"'^'* (2.20) 

It is well known that the canonical action-angle variables are constructed out of the matrix elements 
of T(A,t), and, what is more important for us, this leads to an infinite set of conserved charges. In 
particular, the diagonal element of T(A,t), a{X,t), is independent of time. Then one can obtain all, 
local and non-local, conserved charges as coefficients of asymptotics for /mlna(A): 



Imlna{X) = -^X~''Ak ,|A| — ^00 (2.21) 

k=0 

00 

Imlna{X) = -^X~''Bk ,|A| — ^0 (2.22) 



k= 



where A/^ and iJ^ are global charges given by some known expressions in terms of the action variables 
(see for the details |21H23| ). It is important that these charges can be represented as integrals of 
the /oca/jll i.e., depending on x, quantities: 

Ak = J CLk{x) a-iid Bi^ = J dx bi^{x) (2.23) 

There is a well-established recursive procedure to obtain ak{x) and bk{x), e.g.: 

l S^{x)dSHx)-S\x)dS\x) 

= 4 iTsHx) ^'-''^ 

a^{x) = \riijdS\x)dS\x) (2.25) 
8 

While oq is proportional to the momentum, oi is easily recognized as the Hamiltonian density (j2.7p 
of the model. The higher classical charges could be obtained in the same vein. 

Instead, we will proceed following |24j. As explained in |24) . in the classical theory an arbitrary 
0(1, 2) invarianjf] can be written in terms of three basis invariants and their space derivatives: 

C = S'Si; n{x) = d^S'd^Sf, Q^{x) = e'^''Sid^SjdlSk (2.26) 



*This should not be confused with local charges defined as conserved quantities that locally depend on a field 
S'^ix) and a finite number of its spatial derivatives. In fact, it is well known that series (|2.2ip defines exactly such 
charges, while charges 6^, defined by (|2.22p . have the form of spatial integrals (except for 61 that is just S''^{x) — 1). 

^Actually, in |24| only the 0(3) case was considered. But as we already stressed, both, su{2) and su{l, 1), models 
have exactly the same set of invariants if written in invariant form, which is the case for C, T-L and Q3. 
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where C, T-L{x), and Qsix) correspond to the densities of the Casimir, Hamiltonian and the next 
cubic conserved chargeqj- In Section [3] the appropriate quantum operators wih be derived and, in 
particular, the Casimir operator wih be shown to correspond, after removing regularization, to an 
integrable representation of the su{l, 1) algebra. We will also show that in the quantum theory, 
these fundamental invariants can be written as traces of regularized operators in a simple form. 

2.2 Quantization 

The quantization of the model involves promoting the classical fields to quantum operators 
which satisfy the canonical commutation relations obtained from (12. 3p or (12. Sp . 



[S^ix);S^{y)] =±S^ix)5{x-y) 

(2.27) 

[S-ix)-S+{y)] = 2eS\x)6ix - y) 

and the conjugation conditions: (5'^(x))* = S^{x) and (5^(x))* = S^{x). 

The next step is to construct the representation of the algebra (j2.27p in the ferromagnetic 
vacuum: 

S^{x)\0) = e|0) , 5-(x)|0) = (2.28) 

One possible way is to proceed as in other continuous integrable models U and introduce the 
vectors: 

/n n 
^1 ) ■ ■ ■ ) Xyi )X{S+{x,m (2.29) 
i=i j=i 

where the fn{xi, . . . , Xn) are continuous and sufficiently fast decreasing, symmetric functions of {xn} 
for the integral (j2.29p to be well defined. So, the representation space is given by the linear span 
of the vectors \ fn), n = 1,2, . . .. The scalar product in the n-particle subspace, {gn\fn), is positive 
only in the sn(l,l) case, it is indefinite in the su{2) case |18| . Hence, only in the non-compact 
case it is possible to construct physically meaningful representations of the algebra (I2.27P in the 
ferromagnetic vacuum. Thus, in order to obtain consistent result it is mandatory to set e = 1 in 
(j2. 1112. 281) . We would like to stress that although there is no positively defined scalar product for 
the su{2) case with respect to the ferromagnetic vacuum, it is, however, possible to construct a 
representation with e = — 1 in a non- ferromagnetic vacuum |25| . 
The quantum /^-operator for the isotropic case has the form0 

i ( S\x) -S+{x) 



^(^'^) = aI 5-(x) -S^{x) ) (2.30) 



''For example, the quartic conserved charge Q4 = J dx ^d^S^d'^Si + 5/4 (9a:5"9i^Si)^j can be written as a non-Unear 
function of H and Q3 in the form: Qa = j dx {qH'^ + i^2i£^±«i^ 

'^The £(A, a;)-operator (|2.30p is the quantum version of tlie classical /I'^''^^(A, a;)-operator, which is obtained from 
L^'^^^ (see the footnote on page 6) by the following gauge transformation (we also take A — i- i/a): d(t> = C^'^^\X, x)(f>; 

(j) = Alp; C'-'^''' {X, x) = AL'-'^'^ , where A = ^ *^ ) ' '^^^^ form of /l'-'^'^{\,x) is more convenient in the quantum 

case because the quantum _R-matrix H2.33|l has the standard (as in su{2) case) form. The fact that A is not unitary 
is not a problem, as ip (or (p) does not have any probabilistic interpretation as in the case of Quantum Mechanics. 
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Here A again is the spectral parameter. The monodromy matrix on the interval [x_,rE_|_] is defined 



as: 



r^^+(A) = Pexp / dx C{X,x 



In this case, the bilinear relation, 



(1) (2) (2) (1) 

R{\i - \2)T^^{\i)T^^{\2) = r^_MA2)T^+(Ai)fi(Ai - A2) 

(1) (2) 

where we denote A = A<S) \ and A = \ ® A, holds for the quantum i?-matrix given 



(2.31) 



(2.32) 



(2.33) 



a=0 



with wq{\) = a — y2 and ^1,2,3 = — */2, provided that 



R{\i - A2) 



(1) (2) (1) (2) 

£ (Ai, x) + C (A2, x) + £ (Ai, x) o C (A2, x) 



(1) (2) (2) (1) 

£ (Ai, x) + C (A2, x) + C (A2, x) o C (Ai, x) 



R{\i - A2 



(2.34) 



is satisfied by the quantum £-operator. The o-product was originally introduced in |18| in order to 
make sensible the ill-defined product of operators at the same point, and is defined for any pair of 
local fields A{x) and B{x) as: 



A{x) o B{x) = hm i /""^^ dii di2 A{i^)B{i; 



(2.35) 



Even though (j2.35p is enough for (|2.34p to hold, it is not in general a good prescription for the 
regularized product of operators at the same point, as it takes care of only one type of singularity 
appearing in the (j2.34p . A quick analysis of (|2.35p shows that its action formally corresponds to 
dividing by 6{0) the product of operators at the same point. Namely, it is easy to see that the 
product of two operators satisfying the algebra ()2.27p is regularized by the product ()2.35p in the 
following sense: 



[53(x),5+(x)] =S+{x)6{0) [S^{x)°S+{x)] = 5+(x) 



(2.36) 



Before moving on, we note that only the second expression in (|2.36p is the necessary algebraic 
structure to render (I2.34P true. But now it is clear that (I2.35P is neither suitable for dealing 
with the product of regular operators, nor of operators plagued with more severe singularities. For 
instance, if one takes ^(x) = B{x) = 1, the product defined by (I2.35P implies that lol =0. Hence, 
it is not really a regularization, but a very specific procedure which demands a special care; later 
in section [3l we will show that this problem can be properly addressed by introducing regularized 



^Everywhere below we use the standard su{2) matrices aa = (1,0"^) 
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5-operators. We emphasize that the the vahdity of (|2.32p and, therefore the standard methods of 
the quantum inverse scattering method [ll[26l[27] rehes on the product (|2.35p . 

After passing to the infinite line limit, the operators in the monodromy matrix (|2.3ip satisfy the 
standard commutation relations and the quantities tr [T(A)] form a family of commuting operators, 
which can be regarded as the integrals of motion of the LL-model, and the eigenfunctions of which 
take the form: 

n 

\Xi...\n)=YlB{Xi)\0) (2.37) 

i=l 

with the corresponding eigenvalues given by: 

n 

Y.log[A{X,)] (2.38) 

i=l 

Here, however, one faces an important problem absent in the classical counterpart and the 
lattice version, namely, the difficulty to extract from tr [T(A)] the commuting quantities expressed 
as integrals of local densities, as in (j2. 21112. 2"3]) in the classical case. The reason behind this difficulty 
lies in the fact that the local charges contain operator products at the same point, thus making the 
integrals of motion ill-defined quantities. In particular, one cannot proceed as in the classical case 
and simply decompose the series: 

oo 

tr [r(A)] = J^/„A" (2.39) 

n=0 

to obtain the local charges /„. Hence, the construction of the local integrals of motion poses a 
highly non-trivial problem in the continuous quantum theory, which, from the field theory point of 
view, corresponds to the renormalization procedure (see for example |28|). 

In |18| only the first two charges: the number operator M and the momentum operator V 
were properly constructed. To formulate them, it was necessary to introduce a new set of bosonic 
fields ^nix), corresponding to n-particle clusters, so that they annihilate the ferromagnetic vacuum, 
^n{x)\0) = 0, and satisfy the following algebra: 

= 6mnS{x - y) (2.40) 

In terms of such ^-fields, one can represent the 5-operators as: 

oo 

S'{X) = sl + ^sM{x)^n{x) 

n=l 

oo 

S+ix) = (x) + 4^Ui{x)^n{x) (2.41) 

n=l 
oo 

S~{X) = S^^iix) + S+¥^{x)^n+l{x) 
n=l 

where Sq = 1, Sq = \f2 and s\ = n, = ^n(n + 1) for n > 1. 



^m{x);^\{y) 
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In terms of the ^-fields the number and momentum operator become: 

M = ^n 1 dx^l{x)^n{x) (2.42) 

n=l •' 

oo . „ 

^ = E i / (a,.^t (x)M/„(x) - ^i{x)d^^n{x)) (2.43) 

n=l 

However, the quantum Hamiltonian was not found, only the action of the local quantum- 
mechanical Hamiltonian on one- and two-particle states was guessed. More recently, in [7] it was 
shown that in order to obtain the local quantum Hamiltonian for the arbitrary n-particle sector, it is 
necessary to simultaneously regularize the ill-defined operator product and construct the self-adjoint 
extensions. The first step in this construction is to introduce a regularized quantum Hamiltonian, 
which in [7] was achieved by means of the split-point regularization method. 

Hq = \}^^j dxdyF,{x,y) {-d,S^{x)dyS%y) + d,S+{x)dyS-{y)+ 

+ d,dy [S'^{x)5ix - y)] - d^dy5{x - y)] (2.44) 

where F^{x,y) is an arbitrary smooth and sufficiently fast decreasing symmetric function of (x,y), 
depending on some parameter e, so that the integration in (|2.44p is well-defined and that in the 
limit of vanishing e one has: 

lira Fe{x,y) = 6{x - y) (2.45) 

It is clear from the definition of Hq that it annihilates the ferromagnetic vacuum: 

Hq\0) = (2.46) 

It is also not difficult to compute the action of the Hamiltonian (I2.44P on the general n-particle 
state (f2:29]l 

/n 
[A2 fn{^)]llS+{xm (2.47) 
i=l 

r. n 

+ E / n {mn{^) - 9./n(x)]":5;;: + 9,9,/„(x)u.,=,.J n s'-{x^)\..=.^o) 

where A2 = df stands for the n -dimensional Laplacian. Thus for to be an eigenstate of 

the Hamiltonian Hq, one must require the following matching conditions: 



[9./„,(x) - a,/„(x)]^;i:^:+: + d,d,fn{^%,=., = o , vi > j (2.48) 

Thus, it follows from ^(L^ that 

- A2/„(x) = EnU^) ^ Hq\U) = En\fn) (2.49) 

where En is the energy of the n-particle state. We would also like to point out before moving on 
to the construction of the self-adjoint extensions that, as explained in [7], the matching conditions 
(j2.48p lead to the factorization property of the 5-matrix. 
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For simplicity, here we will only consider the two-particle case and refer the interested reader 
to [7] for the details concerning the construction of the self-adjoint extensions for the general n- 
particle sector. First, one introduces the vector space V generated by vectors of the form 

*-(/(t:)) <-°> 

where h{x) G ^''{R^dx), f2{x,y) E ^^{R''\{x = y},dxdy) and f2{x,y)\-,=y = fi{x) G ^^(R). A 
scalar product in V is defined for any pair £ V as 

{^m) = \ I dxgl{x)h{x)+ f dxdygl{x,y)f{x,y) (2.51) 
The operator H : V ^ V is defined by the following two conditions: 

1. The action of H on f2{x,y) is simply the Laplacian — A2 everywhere in \ {x = y}, i.e., 

H^=( j^f)^^) \ (2.52) 



-^2h{x,y) 

2. H is Hermitian in V with respect to the scalar product (|2.5ip . i.e., 

{H^^^) = {^H^) (2.53) 

These conditions together with the closure of V under H, i.e., H'^ G V , fix the form of h and 
we obtain the hamiltonian action: 

HQ\f2)=( nd.-dy)h{x y)p,t\-dlh{x,x) \ 

which coincides with the formula previously guessed by Sklyanin in |18) and agrees with (j2.48p . 

2.3 A^m-operator hierarchy 

As is well-known, in a classical integrable model the time evolution of a field ^{x, t): dt^{x, t) = 
{Ijyi,^{x,t)}, where an arbitrary conserved charge is chosen as the Hamiltonian, is equivalent 
to the hierarchy of A4m\x, x) operators, satisfying the compatibility conditiorH [3]: 

dtC(''\X,x) = d^Mi^'\X,x) + \Mi^'\X,x),d'^\X,x)\ (2.55) 

where C^^^''{X,x) is the classical limit of the quantum £(A, x)-operator (|2.30p (see the footnote on 
page 8). For example, the standard (L, M)-pair of the classical theory (|2.14p can be obtained from 
the pair (£^^^^(A, x), ^(A, x)). 

For a system defined on the interval [xi,X2] with periodic boundary conditions, the solution has 
the form p6l[29]: 

{trrj^^nA),/:^'^^^(/i,x)} =5,A^^^(a;;A,;u)+ [a^^^(x;A,/.),£^^'^(x,m)] (2.56) 



'We omit everywhere the exphcit dependence on time. 
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where 

(1) / (1) (1) \ 

M^,l{x;X,f,)=tv (rj=HA)r(A-/.)r^^(A)j (2.57) 

and r(A — /i) is the classical r-matrix (see also (|3.49p below) and the operation tr is the trace 
operator in the first space in (E> C^: 

(1) (1) / (1) (2)\ 

tr {A(S) B) = tr { A ■ B \ = {tT A)B 

Now, in the infinite interval hmit the A4m\x,x) operator corresponding to the conserved charges 
Im can be obtained from the series A^i^(a;; A, /x) = fJ-~^ A4m\x, x) . 

In this paper we will be more interested in the quantum operators A^m(A, x), which describe the 
temporal evolution of the quantum transition matrix T^^ (A). For the simpler non-linear Schrodinger 
model the details and explicit relations can be found in |26]. For the LL-model, on the other hand, 
due to the difficulties related to obtaining the local integrals of motion it was only possible obtain 
the quantum 7Wm(A, x)-operators: 

i [/„, £(A, x)] = d^MniX, x)+ : [Mn{X, x), £(A, x)] : (2.58) 

for the first two charges: the number operator Iq = M (j2.42p and the momentum operator Ii = V 
(j2.43p with Mq{X,x) = «o"3/2 and A4i{X,x) = —C{X,x), where : ... : denotes the normal ordering 
with respect to the fields ^nix) [IS]. The relation (|2.58p yields: 

i [/n,T-_+(A)] =: A^„(A,x+)r-_+(A) - T^+{X)Mn{X,x.) : (2.59) 

from which the commutation relations follo^4^: 

[AA,log^(A)] =0 , [M,B{X)] = B{X) (2.60) 
[P,log^(A)]=0 , [V,BiX)] = -jB{X) (2.61) 

In |18) the action of the quantum-mechanical Hamiltonian on the two-particle state, as well as 
the commutation relations (in the infinite interval limit): 

[H,log^(A)] =0 , [n,B{X)] = ^B{X) (2.62) 

leading to the spectrum, were conjectured, so that they reproduce the correct classical limit. In the 
next section we will address this issue in details, and the relations (I2.62P will be properly derived. 
Thus, the importance of the M-operators here is to establish a connection between the inverse 
scattering method, the direct diagonalization, and self-adjointness of the correct local conserved 
charges. 



Here we consider these relations in the isotropic hmit. 
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3 Operator regularization 



In this section we propose a method to deal with all singularities associated with the operator 
product at the same point, and derive a number of results as a consequence. In particular, we will 
obtain the quantum analogue for the three fundamental invariants (|2.26p , derive their spectrum and 
the corresponding M^(\;x,t) operators. We will also explain in which sense this method can be 
considered as an alternative to the discretization procedure. 

As previously discussed, the most challenging problem of the LL-model, as well as all continuous 
integrable models, is the difficulty to express the commuting quantities as integrals of local densities. 
The root of this problem can be traced back to the fact that such local charges contain the product 
of operators at the same point, which in their turn make the aforementioned integrals very ill-defined 
quantities. Our solution to this problem is to deal from the very beginning with a set of regularized 



where J^^{x,y) is some smooth, rapidly decreasing, symmetric function of {x,y), depending on some 
parameter e that renders (13. ip well-defined. More conditions on this J^^{x,y) will be imposed later. 
Our goal will be a complete reformulation of the singular theory (which involves the bare S'-fields) 
in terms of the regularized S'j-(x)-fields. 

It is important to emphasize that the introduction of such a function is more in the spirit of 
the F^{x,y) introduced by [Tj^ to regularize the quantum Hamiltonian (see equation (12.44P ). and 
that it differs in its essence from the approach employed by Sklyanin in |18| . where the product of 
operators at the same point was modified by the introduction of the o-product, defined in (j2.35p . 
Here, the operator product is kept unchanged, i.e., it is still the usual operator product, while the 
operators themselves are regularized according to p.ip . and, as it will become clear below, also 
renormalized. In what follows, we will show that the introduction of such regularization allows us 
to satisfy the bilinear relation (j2.32p without the need to resort to the product ()2.35p . 

As a first step in this direction we obtain the algebra satisfied by the regularized Sjr{x). Inverting 
the relation (|3.ip . 



we introduce G^{x,y) a smooth, rapidly decreasing, symmetric function of {x,y), depending on e 
and such that: 



Now, with the use of the algebra satisfied by the bare operators (j2.27p . one easily derives the 

^^Although the function Ft{x,y) introduced in [7] differs from the function Te{x,y) we use in the definition of 
p.ip . the connection between the two is given, as will become obvious from the text below, by the expression: 
Ft{x,y) ~ J daTt{x,a)Tt{a,y). Although it was not important for the purposes of 0, it will become clear that 
it is necessary to introduce here the T^{x,y) function to deal with singularities in the algebra and the Yang-Baxter 
relation. This will effectively affect the scaling of the operators, and result in correct su{l, 1) representation. 



fields. 






(3.2) 




(3.3) 
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(3.4) 



algebra for the regularized 5jr-operators: 

[S^ixy,S+iy)] =2 J dzlC,ix,y,z)SUz) 
where the kernel }Ce{x,y,z) is the symmetric function of {x,y) defined as follows: 

]C,{x,y,z) = j diF,{x,OJ'e{y,i)Q.{z,i) 

Let us note, that although the algebra of the regularized operators is non-local, it is however local 
after removing the regularization e — > (see the Eq. (j3.10p below). Having obtained the algebra 
for the S'j--operators, we return to the bilinear relation (I2.32p . which we now write in terms of the 
regularized operators S^. As discussed previously, the conditions needed to ensure the validity of 
^(22^ are of the form: 

[S^x)- S^[x)\ = ±S$[x) and [S^{x); S%{x)\ = 2S%{x) (3.6) 

Thus, as long as the algebra ()3.6p is satisfied, the relation (|2.32p follows. It is easy to show that 
this requirement leads to the following constraint on the T^{x,y) function: 



(3.5) 



lim [Jv(x)V 

e-i>0 



lim FAx) 



(3.7) 



We can conclude that with this restriction on the T^[x,y) function in the e — )• limit, the algebra 
(j3.4p reduces to the conditions (|3.6p in the limit x = y, e — )• and, therefore, the relation (|2.34p 
is valid in the e — )• limit, with o taken as the usual operator product, but with the C -operator 
defined in terms of the regularized Sj^- fields, i.e.. 



lim { R{Xi - Aa) 

e-S-O 



(1) (2) (1) (2) 

£jr(Al,x) + Cjr{X2,x) + Cjr{Xi,x) ■ Cjr{\2,x) 



lim ■ 



(1) (2) (2) (1) 

Cjr{\i,x) + £jr(A2,x) + Cjr{\2,x) ■ Cjr{\i,x) 



R{\i - A2) 



and 



Cjr{\x) = - 



SUx) -SU 



'T 

S^ix 



(3.8) 



(3.9) 



-Sjr{x) 

But now, we can use (j3.8p to ensure the validity of the bilinear relation (j2.32p in the limit e — )■ 0, with 
the monodromy matrix (|2.3ip defined with respect to the regularized £jF-operator. As explained 
before, the Eq. (I2.32p is fundamental for the applicability of the quantum inverse scattering method. 
Although the requirement ()3.7p does not have a non-trivial solution in the class of usual functions, 
it does have a solution in the class of generalized functions, and can be written as 

6e{x - y) 



^e{x,y) 



5M 



(3.10) 



where 5^{x) is some regularization of the J-function, so 5e(0) is a well-defined constant before taking 
the limit e — )• 0. Thus, the function J^^{x,y) is not only a regularization, but a renormalization 
as well. Hence, in the case of a product of two operators it will lead to a product like (I2.35p . but 
without its problematic, as previously discussed, behavior on regular fields. 
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3.1 Casimir operator 

One question that naturally arises in this context, is what happens in the quantum theory to 
the classical constraint (12.21): 



C = S^{x)S''^{x)-S+{x)S"{x) = l (3.11) 

It is not difficult to see that the quantum operator defined by C is formally a Casimir operator of 
the algebra (|2.27p . but once again, due to the presence of the product of operators at the same point, 
it is clear that the expression (jS.lip is not well-defined. There are two possible ways to circumvent 
this problem. One can either consider the bare operators and replace the product by the Sklyanin 
o-product (j2.35p . or write ()3.1ip using regularized Sjr operators, but with usual product. Even 
though, we have already pointed out some of the advantages of using the regularized fields when 
compared to the product (j2.35p . we would like to work out this point in more detail for the product 
(I2.35p . so that one of the aforementioned problems of this product becomes clear. 

We already saw (see the discussion after (|2.36p ) that 1 o 1 =0. Here we give an example of an 
even more serious inconsistency of o-product. Noting that: 

rx+A fx+A 



S\x) o S\xm = hm i r dii r dC2 S3(6)5=^(6)|0) 
A->0 A 

1 f-x+A fX+A ^2 

= hm - / d^i / d^2 |0) = lim — 10) = (3.12) 

A Jx Jx A 



one arrives at meaningless result. Hence, we see that working directly with bare operators, even 
with the o-product (|2.35p . does not allow one to perform a completely consistent regularization of 
the model. 

The second possibility is to use the regularized Sjr operators. Let us begin with the general 
expression for a candidate for the Casimir operator: 

= S^x)S^x) - S^x)S^{x) + 7i5|r = K, (3.13) 

Here the parameters 71 and are introduced to reflect the fact that in the quantum theory the 
operators S~^{x) and S~{x) do not commute. We also note, that now both sides depend on the 
regularization parameter e, and, although the expression (I3.13P is not a Casimir of the algebra (13. 4p 
anymore, it can be checked that it becomes one in the limit e — ?• 0. Thus, more care is needed when 
using the quantum regularized constraint ()3.13p . and, in general, one cannot assume that C = 1, as 
in the classical case, since we use rescaled operators. We will now determine the parameter 71, and 
in subsequent sections also the parameter k^. In the e — )• limit it will be shown that we indeed 
constructed spin s = 1/2 representation from the Sjr fields. 

Using the i2j--operator ()3.9p and the algebra (13. 4p . one finds (in the e — )■ limit): 

(£^)2 = 1-1(k,-(7i + 1)5^) (3.14) 

Then, for (Cj^)^ to be uniquely defined, i.e. so that (>Cjf)'^ = (>Cj-) • (Cj^)"^ = (Cj^)^ • (Cj^), one must 
impose the condition [Cjr, Sj] = 0, which leads to the determination of the 71 parameter: 

71 = -1 (3.15) 
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Thus, the regularized constraint (|3.13p takes the form: 

Ct = SUx)SUx) - \ {S%(x)S~^{x) + S^(x)5+(x)] = (3.16) 

We also note the following useful expression, which follows from (|3.14p . and which will be used in 
the subsequent section when deriving the corresponding M operator: 



1 + 4k, . ■ 

^^^^ 



(3.17) 



Finally, we renormalize the ferromagnetic vacuum in the following way. Let us define the ferro- 
magnetic vacuum with the bare operator as before: 

53|0) = ao|0), 5"|0)=0 (3.18) 

which in terms of the S^-p operator takes the form: 

5^10) = C|0), 5^10) = (3.19) 

where the parameters and will be determined at a later stage. Here we only notice, that using 
the definition p.ip . and taking the limit e — )■ 0, the relation between the two parameters should be 
of the form: 

and from the Casimir operator (j3.13p we obtain also the relation between the C,^ and parameters: 

Cl -Ce = K, (3.21) 

3.2 J^-regularized Hamiltonian 

With the regularized Sj^-fields, all the products are now well defined, and we can proceed to the 
construction of the local conserved charges. The first step is to rewrite the quantum Hamiltonian 
(|2.44p in terms of the J-'-regularized fields: 

n = ^ jdx [-2d^sUx)d,sUx) + d,s^x)d^s^{x) + d,s^{x)d^s+{x)] (3.22) 

where 7y(e) is an arbitrary parameter which will be chosen lateil^. We now show that the same 
n-particle states of the form 

/n n 
Y[dxifn{xi,...,Xn)llS+\0) (3.23) 
i=l j=l 



^■^In fact, one could start from a more general form: 



n^^Jdx [~2r,i{e)d^SUx)d.S%{x) + r,2{€)d,S+(x)d^S^{x)+r,s(e)d^S^{x)d^S+{x)] 

where the parameters rii{e), i = 1,2,3 are arbitrary. However, the diagonalization condition will require setting 
771(e) = '72(e) = 773(e) = 77(6), thus, we obtain the Hamiltonian in the form p.22|) . 
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provide a representation space for the su(l, 1) algebra for the operators in terms of Sjr fields. As 
before, the fn{xi, . . . ,Xn) are continuous and sufficiently fast decreasing, symmetric functions of 
{xn} for the integral (|3.23|) to be well-defined. The action of (|3.22|) on the ferromagnetic vacuum 
gives the vacuum energy: 

^|0> = ^^^4^ / dxdzd,dMx,y,z)\,=y\0)=ho{e)\0) (3.24) 



2 

while the action on the one-particle state gives: 



{n - ho{e))\h) ^ I dx dlh{x)S^{x)\0) (3.25) 

Here and below we use the symbol 5{0) for lime_s>o <^e(0). For the more complex two-particle sector, 
it is possible to show that: 

(?^-/io(e))|/2)^-^^ / dxdyA2/2(x,y)5+(x)5+(y)|0) + 



+ ??(e) 



52(0) 

ao 1 



dx [{d. - dy)f2{x,y)]lZ:t':S+{x)S+{x 



<52(0) 25(0) 

dx d^dyf2 (x, y)U=yS+ {x)S+ (x) 1 0) (3.26) 



52(0) 

Thus, choosing the old matching condition (|2.48p we find the value of the ao parameter: 

oo = 1 + (3.27) 

It is easy to understand this result by considering the Casimir operatoi[^ ()3.13p . Remembering the 
relations (I3.20p and (I3.2ip we obtain the parameter in the e — >• limit: 

= SUx)S^j.{x) - \ [S+{x)S^{x) + S-^{x)S%{x)\ = -1/4 (3.28) 

Thus, the constructed sn(l, 1) representation indeed corresponds to spin s = 1/2. This is because 
the S'jr operators were introduced in such a way that at each fixed point x, the su(l, 1) algebra 
they form in the e — )■ limit (I3.6p is well defined. This is unlike the discretization procedure, where 
one needs to introduce discrete 5„ operators and construct non-integrable (i.e., s 7^ Y2,l,ete.) 
representations. 

We can now also fix the value of r/(e), by choosing its limiting behavior: rj{e) — )• Eq^-^ when 
e — )• 0, where Eq is an arbitrary finite constant. Thus, the expressions (|3.25p and (|3.26p take the 
form: 



(n - ho{e))\fi) ^ -Eo j dxdlh{x)S^{x 



{n - ho{e))\f2) ^ -Eo I dx dy A2/2(x, y) S+ (x) S+ {y)\0) (3.29) 



^We emphasize that it is a Casimir only in the e — > limit. 



18 



Consequently, it is possible to reproduce the construction of the self- adjoint extensions for the 
two-particle sector of [7], as discussed in the end of the Section [2.21 More generally, it is possible to 
show that that relations similar to (|2.48p and (|2.49|) hold for the J^- regularized Hamiltonian (|3.22|) . 
so that the general construction of the self-adjoint extension of [7] also goes through in the n-particle 
sector for ([3^2]). 

Finally, we note that the regularized Hamiltonian (|3.22p can now be written in the simple form, 
using the trace and operators as follows: 

n = ho + — dx r/(e)tr [d^C^X, x) ■ d^Cr{\ x)] (3.30) 

This would not have been possible without the introduction of the regularized fields. Thus, our first 
result is that the form of the quantum regularized Hamiltonian (|3.30p essentially coincides with the 
classical fundamental invariant in (12.26p . with classical fields replaced by the S'jr operators. Below 
we will also obtain the other higher-order fundamental invariant in terms of the Sj: fields, and show 
that it has the same form as its classical counterpart. Moreover, it will also be written in the form 
of the trace of £jr-operators. 



3.3 Spectrum 

In the previous section we have found the well-defined quantum Hamiltonian ()3.30p and derived 
its spectrum by means of the direct diagonalization. In this section we show that the spectrum can 
be also obtained from (|3.3Up by means of the inverse scattering method. We establish a connection 
between the two methods by deriving the commutation relations (|2.62p . 

The first step consists of finding the operator A^2(A,x) (see Section ()2.3p ) which satisfies the 
relation: 

i[H;Cr{\x)] = d^M2{\x) + [M2{X,x)-CT{X,x)\ (3.31) 
It is not difficult to compute the left hand side of (j3.3ip : 

i[n,c^{x,x)] ^ --^a,[£^(A,x),a,x^(A,x)] (3.32) 

Now, using the relation (I3.17P it is easy to check that the A^2(A,x) operator has the form: 

M2{\x) = a2Cr{\x) + l32[Cr{\x),d^Cr{\x)] (3.33) 

where the 02 and /32 coefficients have the form: 

I + Ak^ ^ X , . 

"^ = 2A^' ^' = -'m ^^^^ ^ 

Finally, from the equation (|3.3ip one obtains: 

i \%]T^+{X) \ = M2{\x+)T^+{X) -T^, + {X)M2{\x^) (3.35) 
where the transfer matrix is defined with respect to Cjr, and decomposing the Ai2{X,x) as: 

A^2(A, x) = cr+A^+(A, x) + a-M^iX, x) + asMliX, x) (3.36) 
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where: 



MliX,x) = '-^SUx) + § [S+{x)d,Sj 



2/32 



(3.37) 



the equation ()3.3ip is reduced to the set of commutation relations for the elements of the transfer 
matrix (IZSTI) : 

n;Alt{X)] = Mt{X,x+)C-, + {\) + MUx,x+)Alt{X) - Alt{X)MUX,x.) - B'^t{X)M^ 

n;B-,+{X)\ = Mt{X,x+)D-^tiX) + Ml{X,x+)B-^+{X) - A-^t{X)MtiX,x.) + B^^^{X)MU^^ 

(3.38) 

"^;C^'+(A)] =>[2"(A,^+)^^+(A)->(i(A,x+)C,^+(A)-I),^+(A)>[2"(A,x-)-C,^+(A)>[i(A,x_) 
?^;I)^+(A)] =A^2-(A,x+)i?f+(A)-A^i(A,x+)I)^+(A)+I)^-+(A)A^i(A,x_)-C,-+(A)A^+(A,x_) 

Passing to the infinite interval limit is accomplished by introducing the transfer matrix Too (A) 
as the limit: 



Too (A) = lim 

XJ-— ^±oo 



e(-x+,A)r,^_+(A)e(x_,A) 



(3.39) 



where e{x, A) = exp («o'3/ax) . Then, by multiplying both sides of the equations (|3.38p by e{—x^, A) 
from the left and e(x_, A) from the right, taking the limit x± — ?■ zboo and noting that Sjr{x) ^^^^ i 
and S^{x) 



0, one concludes that (in the limit e — t- 0): 




[^;^oo(A)] =0 


(3.40) 


[H;i?oo(A)] = '"f^)i?oo(A) 


(3.41) 


[^;Coo(A)] = -^^^Coo(A) 


(3.42) 


[n;D^{X)]=0 


(3.43) 



Thus, we have proved the relations (I2.62p . from which the spectrum of the quantum Hamiltonian 
follows. This shows how to make a connection between the regularization scheme we have proposed, 
the direct diagonalization, and the inverse scattering method. 



3.4 Higher-order charges 

In this section we construct the quantum version of the cubic in fields invariant in (I2.26P and 
show that, similarly to the Hamiltonian, the density can be expressed in terms of the trace of a 
product of £jr- fields. We then conclude that since all three invariants (I2.26P have such fornJ^ 
and it coincides with their corresponding classical counterpart, any conserved polynomial should be 

^*It is obvious that the Casimir operator can also be written in such form C — J dxtr (Cj^Cj^) 
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expressed in terms of the three fundamental invariants, and, therefore can be written as a product 
of >Cj-- fields and their derivatives. Finally, based on this, we give a general prescription to obtain 
the corresponding A4m operator, which will establish the connection with the inverse scattering 
method. 

Let us first give the final result for the cubic quantum regularized invariant Q3: 

Qs = i^-^^ I dx tr {C^d^C^d^C^) (3.44) 

where the constant co(e) will be fixed later. To derive it, one starts with a general cubic polynomial 
of the form: 

Q3 = icoie) [ dx [aiS'^jrd^S+dlS^ + a2S^jrd^S^dlS+ + asS+d^S^jrd^S^ (3.45) 



3 f>2 0+ 



52 c<3- 



where the coefficients oi, ...,06 are arbitrary and will be fixed from the diagonalization condition. 
It is easy to show, using the algebra (j3.4p . that an arbitrary permutation of the derivatives will lead 
to the same form (|3.45p . Now, after some tedious but straightforward calculations, one can show 
that the vector (I2.29P diagonalizes the Q3 operator (13.45^ . provided two conditions: i) one has to 
fix the relative coefficients oi = — l,a2 = IjOs = 1)04 = — 1)^5 = — l^^e = 1; and ii) to cancel 
the boundary terms, one has to use the same matching condition as before (see (I2.48P ). as well as 
use the oq parameter (j3.27p . which lead to an integrable representation. Finally, with the relative 
coefficients fixed, one can collect the terms into a compact form: 



Q3 = ico(e) j dx e'^\Sr)rd.{Sr)MST)i. 



(3.46) 



From here, one can easily find the result (|3.44p . where we expressed {Sj^)i fields in terms of the 



vCjr(A, j;)-operators. Choosing the function Co(e) such that Co(e) — )• Q, 



4aoC(e) 



when e — )• 0, where Qo 



is an arbitrary finite constant, one can show that, for example, for the two-particle sector: 
Q3I/2) ^-^Qo I dyidy2 {d^y,+dl) f2{yuy2)S+{yi)S+iy2m 



(3.47) 



Collecting the results of the previous sections, namely i) all three quantum fundamental invari- 
ants have the same form as their classical counterparts, when expressed in terms of the regularized 
operators; and ii) the three invariants can be expressed in terms of traces of Cj^{X, x)-operator 
product, we can now claim that an arbitrary conserved charge can be decomposed in terms of these 
invariants and their derivatives. 

We now give a general prescription of obtaining the corresponding TW^-operators, and, thus, 
making a connection to inverse scattering method. Since an arbitrary charge, as we have shown, 
can be expressed in terms of the £jr(A, x)-operator product, it is convenient to use the relation: 



(1) (2) 
£jr(x, Ai),£jr(y, A2) 



dz K^{x,y,z) 



(1) (2) 
r{Xi,\2),£j^{z, Ai) + £t{z, A2) 



(3.48) 
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where the r-matrix has the form: 

r(Ai,A2) 



2(Ai - As) 



(3.49) 



Note, that in ah final expressions, in the Ai — t- A2 = A hmit, the singularity V('**i~'^2) will cancel. 
Then, for example, for the Hamiltonian, one finds: 



[?^,£^(y,A2)] =tr (i) J 



dx 



(1) (2) 

^(X, Ai),£jr(?/, A2) 



which in the limit e — )■ and Ai — )• A2 gives: 



Ar?(€) 
25M 



dy [Cjr, dyCjr] 



(3.50) 



(3.51) 



leading to the previous result (13 .3111313411 . 

Similarly, using (j3.48p . one can derive the corresponding Mm operator for any charge, even 
though it becomes very tedious with each higher order. For example, for Q3, a similar calculation 
results in the following expression: 



(3.52) 



where A{Cj^) is some function of the operator Cj^{y, A). 



4 Self-adjoint operators and extensions 

The final check of our construction is the verification of the self-adjointness of each operator. 
The self-adjointess of the Hamiltonian, together with the construction of self-adjoint extensions, 
was first worked out in [7] (see also Section fl2.2p ). Here we will also generalize it to the next order 
fundamental invariant Q3. 

Consider the scalar product (see [7]): 

{~9n\fn) = E e-~''-Cp J d^^-C(x)/„(x)||,^,^„}=,„ (4.1) 

Partitions 

where {X^} is a partition of the set {xi}: 

Mp 

1^ Xm = {X'i} and Xm Xn = SmnXm 
m=l 

tm is a 'collective' coordinate for all Xi G Xm and Cp are positive combinatorial factors. We note 
that the case with n = 2 is equivalent to the scalar product (j2.5ip . We consider here in details the 
simpler two-particle sector with the scalar product ()2.5ip . and in the end comment on the general 
n-particle case. Introducing the action of Q3 on I/2) as: 
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with A3 = u) XlILi ^1^1 where n denotes the appropriate dimension, which in this case is n = 2 and 
uj G iM. Self-adjointness with respect to (j2.5ip . {Q-ig2\fn) = {g2\Q2,fn) , demands 



dxgl{x)hfi{x) - / dx dy g2{x,y) A3 f2{x,y) 



\ j dx [i^giix]^ ^^^^^ ~ ^^'^^ i^3g2{x,y))* f2{x,y) (4.3) 



2 

Then, we compute 
/ dxdygl{x,y)/\2,f2{x,y)= j dx dy {A3g2{x,y))* f2{x,y) + 

Jx^y Jx^y 

+ m.,v) (31 - Si) U..y) + (8? - 8})m^.y)U^.y) - 

- a,r^(x. y)d,h{x, y) + d,g',(x, y)d,k{x. (44) 

which is equivalent to 



3 



dx {g2{^^y)[^{9l-dy)] f2ix,y)}UJ^l- / dx dy gl{x,y)A3f2{x,y) 

x+y 



\b}{dl-d'^^~g2{x,y)X h{x,y)^ \ dx dy {A3g2{x,y))* f2{x,y) (4.5) 

y—^ ^ J x^y 

Comparison with (|4.3p yields: 

hhix) = 3 [co {dl - d'y) f2{x,y)]lZT_l + i'sMx) (4.6) 

where is any self-adjoint operator in (M, dx). Next, the condition that the space generated by 
1/2) is closed under the action of Q3 allows us to fix the form of ig. 

-^3f2ix,y)\x=y = hh 

= 3[a; {dl-dl)f2{x,y)]'=J^^_\ + l'3h{x) 

-udlf2{x,x) = 3oj (a, + dy) [(9, - dy) f2{x,y)\lz:tl - d,dyf2{x,y)\\,=y] + i'^Mx) (4.7) 

which upon the use of the matching condition (|2.48p for n = 2: 

[dxf2{x,y) - dyf2{x,y)fyZltl = da:dyf2{x,y)\x=y 

leads to the conclusion: 

i'3Mx) = -ojdlf2{x,x) (4.8) 
and finally, we are allowed to conclude that 



0\f)=( ^^m-d'y)f2{x,y)]lzT_l-^d'xf2ix,x) 



(4.9) 



is a self-adjoint operator. It is important to stress that (I4.9p does not impose any additional 
condition on the Hilbert space but rather is a consequence of the only matching conditions ()2.48p 
defining our Hilbert space. This result can be generalized for an arbitrary n-particle sector using 
the S-matrix factorization property (see [7] for details). 
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5 Conclusion 



In this paper we have proposed a method to quantize continuous integrable models on the 
concrete example of Landau-Lifshitz model, without resorting to discretization schemes. We also 
explained how to obtain the quantum trace identities, and established a connection with the inverse 
scattering method. The reason behind quantizing the system directly in the continuous framework 
is that discretization is not always an obvious procedure, and for more complicated models it is very 
desirable to learn to quantize integrable systems directly. The difficulty in quantizing the continuos 
model directly is the not well defined operator products, thus, leading to operator ordering problem. 
Even for the simplest models, such as the non-linear Schrodinger model, the normal ordering is 
a non-trivial procedure in quantum trace identities, and generally relies on some discretization 
procedure. 

For the LL model the situation is more involved, since it is nor obvious how to sort three fields S"^, 
and S~ . One can try to resolve the constraint and use the normal ordering for the unconstrained 
fields, but this approach suffers from various problems, i.e., inability to find the quantum charges, 
etc. Besides, it is very desirable to obtain the quantum charges in terms of the original fields, 
namely, to obtain the quantum analogues of the classical expressions. In the theory with more than 
three fields, even after resolving the constraint, if it is even feasible, one will still face the problem of 
operator ordering. In this work we made the first steps to address all these problems. Although we 
work out these ideas on the particular LL model, and there are still several questions that need to 
be understood and generalized, we believe our method can be useful for quantizing other continuous 
integrable models. 

Our method relies on several seemingly different aspects of the integrability, which we have 
shown to be inter-connected in the end. Firstly, the direct diagonalization requires two different 
procedures: i) operator regularization, and ii) construction of the extended Hilbert space. The 
operator regularization we have introduced is very general, and only satisfies a self-consistency 
condition following from the intertwining relation. This regularization, which at the same time was 
shown to renormalize the operators, has significant consequences in various aspects. Firstly, the 
algebra of the quantum operators was shown to be non-local even for the isotropic model. Let us 
remind, that previously such non-local algebra (Sklyanin algebra) appeared only in the anisotropic 
model. We also emphasize that our regularization is essentially different from the one originally 
introduced in |18H20t [30]. While in the original version the regularization (Sklyanin product) was 
chosen to satisfy the intertwining relation, it regularizes only one type of singularity. In the version 
we propose, the J^-regularization regularizes all singularities at the same time, and leads to well- 
defined operator products, as well as a well-defined algebra. We were able to show that this leads 
to integrable representations, unlike the discretization procedure, where the process of passing from 
discrete to continuous operators requires construction of non-integrable representations. 

Secondly, the singular nature of the LL model (which is why one has to introduce operator 
regularization), requires the construction of self-adjoint extensions. While normally this is being 
ignored in physical models, it is nevertheless crucial for this model, as well as a number of other 
interesting models, and is relevant for strings in particular. To name a few simplest ones, the 
principal chiral model in the RxS^ subsector, the higher-order corrections (see |31)). or the fermionic 
AAF model will as well require the construction of self-adjoint extensions. One of the problems 
we would like to undertake is the nature of the minimal extension to accommodate the entire 
superstring on AdS^ x S^. This is not an obvious question and will require a more detailed study 
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of applications of functional analysis on superspace. 

Together, the operator regularization, subject to quantum integrability, and the construction 
of self-adjoint extensions lead to several consequences: the S-matrix factorization, the nature of 
the representation, the connection to the inverse scattering method, etc. Let us once more stress 
on the latter point: the relation to inverse scattering method is established via the hierarchy of 
Mm operators. Namely, only after introducing the regularized operators and verifying the self- 
adjointness of the operators, it is possible to derive strictly the corresponding Mm operator. We 
have shown this directly on the examples of the quantum Hamiltonian and the higher-order cubic 
invariant. It was not possible to derive the M2 operator corresponding to the Hamiltonian before, 
while in our method, deriving the corresponding Mm operator for any conserved charge, although 
very tedious, is nevertheless very straightforward. 

Finally, we have obtained the three fundamental invariants in the quantum theory such that, 
as in the classical theory, any conserved charge can be decomposed into these three invariants. 
Moreover, we have shown that all of them can be written in terms of the traces of the Cjr operator 
products and their space derivatives. This is exactly how one can obtain the Mm operators. 

There are several problems one can undertake in the future. Even for the LL model, we have 
considered only the isotropic case. We have emphasized throughout the paper that anisotropy 
introduces not only technical difficulties, but is essentially different for the integrability. Namely, 
to satisfy the intertwining relations, one needs to introduce an additional quantum operator, thus, 
modifying the algebra of observables. This gives rise to quadratic Sklyanin algebras, which in 
turn involves the Sklyanin product extensively discussed in the text. It will be very interesting 
to investigate the anisotropic case using our prescription and shed some light on the origin of such 
quadratic algebras and their non-local nature. We expect that the anisotropy will completely change 
the extended Hilbert space and require more intense study of self-adjointness of the operators. We 
only mention the importance of the Sklyanin algebra, as the discretization scheme, which leads 
to the /ocal charges in the discrete version, depends on representations of such algebras (see for 
example [2l[3l[5]). We hope to report on progress in this direction soon. 
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